On the finite section method for computing exponentials of doubly-infinite skew-Hermitian matrices
نویسندگان
چکیده
Computing the exponential of large-scale skew-Hermitian matrices or parts thereof is frequently required in applications. In this work, we consider the task of extracting finite diagonal blocks from a doubly-infinite skew-Hermitian matrix. These matrices usually have unbounded entries which impede the application of many classical techniques from approximation theory. We analyze the decay property of matrix exponentials for several classes of banded skew-Hermitian matrices. Then finite section methods based on the decay property are presented. We use several examples to demonstrate the effectiveness of these methods.
منابع مشابه
Computing the Matrix Geometric Mean of Two HPD Matrices: A Stable Iterative Method
A new iteration scheme for computing the sign of a matrix which has no pure imaginary eigenvalues is presented. Then, by applying a well-known identity in matrix functions theory, an algorithm for computing the geometric mean of two Hermitian positive definite matrices is constructed. Moreover, another efficient algorithm for this purpose is derived free from the computation of principal matrix...
متن کاملCanonical forms for doubly structured matrices and pencils
In this paper we derive canonical forms under structure preserving equivalence transformations for matrices and matrix pencils that have a multiple structure, which is either an H-selfadjoint or H-skew-adjoint structure, where the matrix H is a complex nonsingular Hermitian or skew-Hermitian matrix. Matrices and pencils of such multiple structures arise for example in quantum chemistry in Hartr...
متن کاملAn iterative method for the Hermitian-generalized Hamiltonian solutions to the inverse problem AX=B with a submatrix constraint
In this paper, an iterative method is proposed for solving the matrix inverse problem $AX=B$ for Hermitian-generalized Hamiltonian matrices with a submatrix constraint. By this iterative method, for any initial matrix $A_0$, a solution $A^*$ can be obtained in finite iteration steps in the absence of roundoff errors, and the solution with least norm can be obtained by choosing a special kind of...
متن کاملQuantitative Estimates for the Finite Section Method
The finite section method is a classical scheme to approximate the solution of an infinite system of linear equations. We present quantitative estimates for the rate of the convergence of the finite section method on weighted `-spaces. Our approach uses recent results from the theory of Banach algebras of matrices with off-diagonal decay. Furthermore, we demonstrate that Banach algebra theory p...
متن کاملOptimal Parameter in Hermitian and Skew-Hermitian Splitting Method for Certain Two-by-Two Block Matrices
The optimal parameter of the Hermitian/skew-Hermitian splitting (HSS) iteration method for a real 2-by-2 linear system is obtained. The result is used to determine the optimal parameters for linear systems associated with certain 2-by-2 block matrices, and to estimate the optimal parameters of the HSS iteration method for linear systems with n-by-n real coefficient matrices. Numerical examples ...
متن کامل